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Abstract 

In this paper, optimum trajectories in Earth Transfer Orbit (ETO) for a lunar transportation system are 
proposed. This paper aims at improving the payload ratio of the reusable orbital transfer vehicle (OTV), 
which transports the payload from Low Earth Orbit (LEO) to Lunar Low Orbit (LLO) and returns to LEO. In 
ETO, we discuss ballistic flight using chemical propulsion, multi-impulse flight using electrical propulsion, 
and aero-assisted flight using aero-brake. The feasibility of the OTV is considered. 


1. Nomenclature 

a = semimajor axis [km] 

a dr = radius direction component of thrust [km/s 2 ] 

a de = argument direction component 

of thrust [km/s 2 ] 

a dh = normal direction of the orbit plane 

component of thrust [km/s 2 ] 
a T = thrust acceleration magnitude [km/s 2 ] 

e = eccentricity [-] 

/ = true anomaly [rad] 

h = altitude [km] 

I = specific impulse [sec] 

/ = mean longitude [rad] 

M = mean anomaly [rad] 

n = mean motion [rad/s] 

p = semi-latus rectum [km] 

r = orbital radius [km] 

r a = altitude of apogee [km] 

r p = altitude of perigee [km] 

Q in = heating rate on stagnation point [W/m 2 ] 

Q out = radiation energy [W/m 2 ] 

R e = radius of the Earth [km] 

R n = radius of curvature of nose [m] 

T max = temperature of stagnation point [K] 

t f = total flight time [sec] 


t p = orbital period [sec] 

w = argument of periapsis [rad] 

s = emissivity [-] 

(p = in-plane thrust steering angle [rad] 
p = density of atmosphere [kg/m 3 ] 

cr = coefficient of Stefan-Boltzmann 

[W/(m 2 K 4 )] 

2. Introduction 

The moon is recognized as an important destination 
for space science and exploration. To explore the 
moon efficiently, it is conceived that the special 
spacecraft is used in each segment such as LEO, 
Lunar/Earth Transfer Orbit (LTO/ETO), and LLO. 
For example, the OTV is the transportation 
spacecraft used in LTO/ETO. This paper aims at 
improving the payload ratio of the reusable OTV. 
The reusable OTV system is roughly composed of 
the payload module, the propulsion module, and the 
fuel cartridge. The payload module (included man) is 
separated from the OTV in LTO and is reached LLO 
by its single injection. The other modules (not 
included man) are returned to LEO through ETO. 

In this research, we focus on the optimization 
problem of return trajectories with small fuel 
consumption, because the return flight in ETO is not 
constrained the flight time and the safety compared 
to the flight in LTO. In ETO, we discuss ballistic 
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flight using chemical propulsion multi-impulse flight 
using electrical propulsion, and aero-assisted flight 
using aero-brake. The feasibility of the OTV is 
considered. 



Fig.l Image of Orbital Transfer Vehicle 

3. Lunar Free Return Trajectory 

A lunar free return trajectory is used a return reference 
orbit in ETO. The lunar free return trajectories have 
been studied and used since the Apollo era[l,2]. The 
trajectories have many beneficial characteristics for the 
OTV concept. For example, one characteristic is that 
the spacecraft is automatically returned to the Earth by 
only a single injection maneuver at the Earth departure 
and no further large delta-V is required except for small 
adjustment maneuver thereafter. Thus, the trajectories 
are very easy to operating the system such as the 
payload separation and very useful for saving the fuel 
consumption included the total operation. 

The reference orbit in ETO is calculated using the 
patched-conic approximation. 



Lunar Sphere of Influence 


Fig.2 Lunar Free Return Trajectory 


4. Optimum Trajectories using Low-Thrust in 
ETO 

The optimum trajectories are designed from the 
periapsis in the elliptic orbit after the lunar flyby (Fig. 
3). If the flight time is not constrained, the optimum 


control on the minimum fuel consumption is only the 
continuous injection near the periapsis. It is difficult to 
apply the optimum control the OTV system, because 
the flight time is very long. So, it is assumed that the 
OTV in ETO have the fuel to an extent and can always 
keep the thrust constant. It is also assumed that the 
OTV mass is constant, because the fuel consumption 
caused by low-thrust propulsion is very small compared 
to the initial mass. Then, the optimization problem on 
the minimum fuel consumption is equal to the 
optimization problem on the minimum flight time. 


Earth 


Moon 



Fig.3 Section for Optimum Trajectory Design 
4.1 Configuration 

Now, the optimum control problem is to find the 
control history, u(t), which minimizes the following 
objective function J. Firstly, we consider the two-body 
problem in the two dimensions(in-plane) to simplify the 
problem. 

J = f 0 dt (1) 

subject to the equation of motion in polar coordinate 
system (Gauss’s form of the Lagrange’s planetary 
equations)[2] 

da 2a 2 f p 'j 

- = —\es lnf a dr +-a de J 

< ^ = |{ psinfa dr +[(p+r)cosf + re]a M } (2) 

dl ae r r r \ . r i 2br 

p = a( 1 - e 2 \ r = — , h = nab 

{ ' > 1 + ecosf J 

where is the mean longitude defined the following. 

l = Q + w + M (3) 


In the case of the in-plane motion, / = m • 
the control constraints 
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[ a dr = a T sin(p 
\a d Q = a T cos(p 


(a T = 10~ 7 km/s^j 


( 4 ) 


where the steering angle q> is defined in Fig.4. 


the boundary conditions 

a(0) = a o = 273024 

< e(0) = e 0 = 0.9759 (5) 

l(0)=lg=0 
K*85R e ' 

W=R e + 200 y 

a{tf} = a, = 6778 

< e ( tf )= e ,=0 ( 6 ) 

l(tf) = l t = 2 k 

[\ =r p, r =R e +200 ) 



The trajectory optimizes in the orbital period, and 
each of the optimum trajectories in the period is 
connected. Then, it is considered that the connected 
trajectory is the optimum trajectory in ETO. 


based on calculus of variation, dealing with the 
variables directly, this method can be easily applied to 
various optimum control problems. Especially for the 
problems with state inequality constraints this is 
probably the most powerful numerical solver, since the 
method based on calculus of variation needs much 
effort to solve this kind of problem. 

On the other hand, this method needs a lot of memory 
and long calculation time. But in this point the 
remarkable development of computers in these days has 
overcome these difficulties. The nonlinear problem is 
solved by Sequential Quadratic Programming(SQP) 
method. 

4.3 Numerical Simulation in Quasi Circular Orbit 

Firstly, we discuss the quasi circular orbit near the 
target circular orbit. The transfer time is on the order of 
hundreds of days, which results in thousands of orbital 
revolutions. So, it is difficult to optimize the trajectory 
in total flight time simultaneously. The trajectory 
should be optimized in the orbital period using DCNLP. 
The numerical condition is the following, 

the objective function 

J = fo'dt (7) 

the equation of motion (2) 

the control constraints (4) 

the boundary conditions 


a (0) = a 0 = (R e + 200 ) + ^ ~ ^ = 6678 

< e(0) = e 0 = Va ° ~ rp ° = 0.014975 (8) 

2a 0 

i(o) = i 0 = o 

( r a, - r p. =200 ) 


4.2 Method of Optimization (DCNLP) 

There are a lot of numerical methods for trajectory 
optimization [3]. In this paper, we approximate state 
variables by cubic polynomials, interpolate control 
variables linearly and use collocation to satisfy the 
differential equations. This method is called Direct 
Collocation with Nonlinear Programming (DCNLP). 

This is method which converts the original optimum 
control problem into Nonlinear Programming (NLP) by 
discretizing the originally time continuous functions, 
and solves this NLP [4, 5]. Compared to the method 


a { t p) = a t, = f ree 
< e { t P ) = e t, = f ree 
l(tp) = l t = 2 k 

These equations are normalized by using proper 
characteristic length L, characteristic time T, and so on. 
The results are shown in Fig. 5 and Fig. 6. 
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Fig.5 Optimum trajectory in quasi circular orbit 
(Inertial coordinate system) 
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Fig.6 Optimum thrust steering angle profile 


a(0) = a 0 =(R e + 200)+ r ^p^ = 16578 

< e(0) = e 0 = r ^ r ^ =0.60321 (10) 

2a 0 

i(o)=i 0 =o 

(r a -r p =20000) 

a { t p) = a t, = f ree 

<e (tp) = e , = fr ee (9) 

l(tp) = l t = 2ft 


These equations are normalized by using proper 
characteristic length L, characteristic time T, and so on. 
The results are shown in Fig. 7 and Fig. 8. 
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Fig.7 Optimum trajectory in elliptic circular orbit 
(Inertial coordinate system) 


4.4 Numerical Simulation in Elliptic Orbit 

Next, we discuss the elliptic orbit. The trajectory 
should be also optimized in the orbital period using 
DCNLP. The numerical condition is the following, 

the objective function (7) 

the equation of motion (2) 

the control constraints (4) 

the boundary conditions 
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Fig.8 Optimum thrust steering angle profile 


Fig.9 Optimum steering law for minimum a 


4.5 Analytical Optimum Steering Laws 
(In-plane motion) 

In two dimensions, the semimajor axis and the 
eccentricity must be changed to reach the target LEO. 
To minimize the total flight time for the orbital transfer 
and meet the boundary conditions of the desired final 
orbit, the optimum thrust direction time history must be 
computed. 

Firstly, the control law to minimize the time rate of 
change of the semimajor axis is considered, and is 
derived from the governing equation of motion for a . 

— = \ esinfa T sincp + — a T cos(p ] (1 1) 

dt h \ r J 


The optimum in-plane steering angle that maximizes 
the time rate of change of a is found by steering the 
partial derivative da/ dtp equal to zero. 


!r~ = ~^{^ es i n f a T cos( p - ~ a r s ^ n( p j = 0 

Therefore, the steering law for minimum a is 

esinf 


(p a = tan J 


ecosf + 1 


( 12 ) 


( 13 ) 



Fig. 10 Optimum steering law for minimum d 
(Inertial coordinate system) 


Next, the control law to minimize the time rate of 
change of the eccentricity is considered, and is derived 
from the governing equation of motion for e . 


de 

dt 


- ( psinfa T sin(p + \{p + r) cosf + re] a T cos(p ) 


(14) 


Again, the partial derivative is set to zero. Therefore, 
the steering law for minimum e is 


(p e = tan J 


psinf 


(p + r)cosf + re 


( 15 ) 
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Fig. 11 Optimum steering law for minimum e 

(, a = R e +200 ) 



Fig.13 Optimum steering law for minimum e 
(Inertial coordinate system) 

(a =16578, e = 0.60321) 



Fig.12 Optimum steering law for minimum e 

( e = 0.5) 


It is found that the steering law for minimum e 
depends on the eccentricity and doesn’t depend on the 
semimajor axis from Fig. 11 and Fig.12. 

Finally, the control law to minimize the time rate of 
change of the height of apogee is considered, and is 
derived from the governing equation of motion for r a . 


* i 
dt 


(l+e)— + a 
K J dt 


de 

dt 


(16) 


Again, the partial derivative is set to zero. Therefore, 
the steering law for minimum f a is 


(p r = tan 


rsinf \2ae(l + e) + /?} 
r(p + r)cosf + 2a(l + e)p + r 2 e 


(17) 
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Fig.14 Optimum steering law for minimum r a 

(a = R e +200) 



Fig. 15 Optimum steering law for minimum r a 

(e = 0.5) 

It is found that the steering law for minimum r a also 

depends on the eccentricity and doesn’t depend on the 
semimajor axis from Fig.14 and Fig. 15. 

To determine if each equation (Eq.(13), (15), and 
(17)) minimizes or maximizes the value (a,e, or f a ) 

is implicitly used the following boundary conditions. 

<P a (f = 0) = cp e (f = 0) = (p ra (f = 0) = n (18) 


4.6 Analytical Optimum Steering Laws 
(Out-of-plane motion) 

The control law to minimize the time of rate of change 
of the inclination is considered, and is derived from the 
governing equation of motion for i . 

di _ rcos(w + f) 

It h — a “ h (18) 

The partial derivative is set to zero. Therefore, the 
steering law for minimum i is 

W 2 (if cos(w+f)>6) 

<Pi = u / , x ( 19 ) 

{ 37 /2 (if cos(w+f)<0) 

When the angle w + f is near the ±90 deg, any 
out-of-plane steering is essentially wasted. 

5. System Design for OTV 

This paper aims at improving the payload ratio of the 
reusable OTV. The optimization problem is the 
mini mization of the fuel consumption. In the case of 
the constant thrust, the optimization problem is equal to 
the minimization of the total flight time. To minimize 
the end time for orbital transfer, the optimum thrust 
direction time history must be considered. The 
optimum thrust direction for the minimization of the 
total flight time is investigated. In section 4.5, the three 
optimum steering laws for each parameter are obtained. 

In our approach, each steering law is selected in the 
stage of the return trajectory for the OTV, because the 
operation of the OTV system should be simple. 

Firstly, the initial value problem is computed for 60 
days by using each optimum steering law, and the each 
result is compared. The state equations are numerically 
integrated by using a standard fixed step (lsec), 
fourth-order Runge-Kutta routine. The results are 
shown in Table.l and Fig. 16-19. 


Table.l Comparison among results of 

each steering law (after 60 days) 


Steering law 

a[km] 

e 

rjkm] 

r p [km\ 

<Pa 

35760 

0.8227 

65180 

6340 

Ve 

35860 

0.8204 

65280 

6440 

<Pra 

35770 

0.8210 

65140 

6410 
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Fig. 16 Trajectory using steering law for 60 days 



Fig. 17 Change of eccentricity using steering law 
for 60 days 


Fig.18 Change of apogee and semimajor axis 
using steering law for 60 days 



Fig.19 Change of perigee using steering law 
for 60 days 

Then, in the case of the constant altitude of the 
perigee, the minimization of total flight time is equal to 
the minimization of the altitude of the apogee. Now, the 
altitude of the perigee is not constant, and it is difficult 
to decide that which steering law is the best. However, 
when the altitude of the perigee is nearly approximated 
constant, the optimum steering law is the law for 
minimum a or for minimum r . There is not much 

a 

difference between the minimum a law and the 
minimum r a laws, but the optimum steering law in the 
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first stage is selected the minimum a law, because it 
is easy to operate from the Fig. 10 in Section 4.5. 
However, the change of the perigee should be 
constrained, because the heating rate of the OTV 
caused by atmosphere is increased. It is assumed that 
the OTV system can’t control the lift/drag ratio. Then, 
the operation in the next stage is to raise the altitude of 
the perigee. Now, when the altitude of the perigee is 
about 100- 15 0km, the operation for the perigee up is 
done. For example, when the altitude of the perigee is 
smaller than 100 km, the thrust acceleration direction is 
controlled the velocity direction in the vicinity of the 
apogee. 

Next, when the eccentricity or the semimajor axis is 
small, the optimum steering law is selected the 
minimum r a law, because optimum steering angle 

profiles in the numerical simulation results of the case 
is similar to the minimum r a law from the Fig. 6, Fig.8, 

Fig. 14, Fig. 15. 

Finally, when the semimajor axis is equal to the target 
orbit, the optimum steering law is selected the 
minimum e law. 

We compare the low-thrust using the optimum thrust 
control strategy with the ballistic flight using the 
chemical propulsion. 

Table.2 Comparison low-thrust with 
ballistic flight for OTV 


Flight 

t f [days] 

A m[kg](m i =5000 kg) 

Ballistic 

- 

3250 

Low-thrust 

(optimum) 

410 

1070 


When the target inclination is different from the initial 
inclination, the operation to change the inclination must 
be done. It is considered that the operation is done after 
the arrival of the target orbit (LEO) or these steering 
laws are blended so that the orbital elements are 
controlled simultaneously, In our approach, the 
operation is done in the vicinity of the apogee in the 
first stage, because it is simple and easy to do the next 
operation. 

6. Aero-assist (Atmospheric drag) 

The effect of the atmospheric drag should be 
considered, because the altitude of the perigee is low. 
Firstly, the aero-brake is considered. When the lift/drag 
ratio can be changed, some optimum trajectories are 


considered [7,8,9]. However, it is assumed that the 
OTV system can’t control the lift/drag ratio and the air 
drag acts only in a direction opposite to the velocity 
vector, because it is simple and easy to operate. It is 

assumed that the ballistic coefficient is 
constant 0.02m 2 /kg. The acceleration of the 
air drag is nearly equal to the electric 
propulsion in the altitude of 200km. 

Next, the effect of the aerodynamic heating is 
considered. The top of the OTV is generally reached 
the maximum temperature. The heating rate on the 
stagnation point is calculated using the following 
equation[10]. 

1.825 xl0~ 4 rj7\ T/ 3 / 0rVk 

Qm = rr <p\ h y ( 2 °) 

\J k N 

The loss of the thermal energy is given by the 
following equation called the Stefan-Boltzmann law 
of radiation. 




(a = 5.67* It r*) 


( 21 ) 


The T max is calculated using Q jn = Q ouj . The velocity 
of the OTV is the maximum at the perigee. 


f 1.825*1 O' 4 


/ 


pv m 


( 22 ) 


The model used for air density is a 1976 Standard 
Atmosphere. It is assumed that the value of emmisivity 
is 0.65 and the radius of curvature of the nose is 2m. 
The maximum velocity in the return trajectory is 
10.9406 km/s at the initial perigee. Then, the maximum 
temperature of the OTV is 136 IK at a height of 100km. 
It is possible to make a flight at a height of about 
100km, because the allowable temperature limit of 
tempered glass is about 1800K. 

7. Conclusion 

The optimum return trajectory of the OTV is proposed. 
In this research, it is assumed that it is not included man 
in the reusable modules of the OTV. Then, we focus on 
the optimization problem with small fuel consumption. 
The control strategy is mainly divided into the three 
stages. The steering law is selected the minimum a 
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law in the first stage, the minimum r a law in the 

second stage, and the minimum e law in the final 
stage. When the perigee must be raised or the 
inclination must be changed, the steering law is 
selected the perigee-up law or the minimum i law. 
The problem, which is when it should be switched the 
minimum a law to the minimum r a law, is left. It 

this investigation, it is switched the minimum a law 
to the minimum r law in the altitude of about 50,000 

a 

km. 

The feasibility of the reusable OTV with solar electric 
propulsion is considered. The fuel consumption is 
greatly decreased using the low-thrust (Hall thruster). 
The weight of the OTV is increased by equipping solar 
electric propulsion with, however the increase of mass 
is a few hundreds kg. So, the improvement of the 
payload ratio is the remarkable. The flight time is also 
acceptable. 
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